A Measurement of Primordial Non-Gaussianity Using WMAP 5- Year Temperature 

Skewness Power Spectrum 
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We constrain the primordial non-Gaussianity parameter of the local model /nl using the skewness 
power spectrum associated with the two-to-one cumulant correlator of cosmic microwave background 
temperature anisotropies. This bispectrum-related power spectrum was constructed after weighting 
the temperature map with the appropriate window functions to form an estimator that probes the 
multipolar dependence of the underlying bispectrum associated with the primordial non-Gaussianity. 
We also estimate a separate skewness power spectrum sensitive more strongly to unresolved point 
sources. When compared to previous attempts at measuring the primordial non-Gaussianity with 
WMAP data, our estimators have the main advantage that we do not collapse information to a single 
number. When model fitting the two-to-one skewness power spectrum we make use of bispectra 
generated by the primordial non-Gaussianity, radio point sources, and lensing-secondary correlation. 
We analyze Q, V and W-band WMAP 5-year data using the KQ75 mask out to i max = 600. Using 
V and W-band data and marginalizing over model parameters related to point sources and lensing- 
secondary bispectrum, our overall and preferred constraint on /nl is 11.0±23.7 at the 68% confidence 
level (—36.4 < /nl < 58.4 at 95% confidence). We find no evidence for a non-zero value of /nl even 
marginally at the la level. 

PACS numbers: 98.70.Vc, 98.80.-k, 98.80.Bp, 98.80.Es 



I. INTRODUCTION 



The inflationary paradigm has deservedly become a 
cornerstone of modern cosmology [l|, H, H, 0, H, @| . Infla- 
tion solves the flatness, horizon and the monopole prob- 
lems of the standard Big-Bang cosmology. Furthermore, 
inflation is the prevailing paradigm related to the origin 
of density perturbations that gave rise to the large-scale 
structure we see today. It posits that a nearly exponen- 
tial expansion stretched space in the first moments of 
the early universe and promoted microscopic quantum 
fluctuations to perturbations on cosmological scales to- 
day 0, H, HI OH El • Inflation makes detailed predictions 
for key statistical features of these fluctuations. These 
predictions have now begun to be tested by a range of 
cosmological observations, including cosmic microwave 
background (CMB) temperature anisotropy and polar- 
ization. 

Recent measurements of the CMB with a variety of 
ground, sub-orbital, and space-based experiments have 
provided some of the most stringent tests of inflation 
(e.g., Specifically among the generic predictions 

of inflation, recent CMB measurements with the temper- 
ature anisotropy power spectrum and polarization have 
established (1) a nearly flat geometry, (2) a nearly scale- 
invariant spectrum at large angular scales, (3) adiabatic 



fluctuations, and (4) super-horizon flucuations through 
the temperature-polarization cross spectrum. One major 
prediction of inflation yet to be verified is the stochastic 
background of primordial gravitational waves [pi , [Til. [To] ] . 
While strong limits are expected from Planck [lq |. a de- 
tection of the gravitational wave background is the main 
focus of a next-generation space-based CMB experiment 
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Some other tests of inflation involve the probability 
distribution function and isotropy of the density pertur- 
bations generated by inflation. In the standard slow-roll 
inflationary model the inflaton, the hypothesized scalar 
field or particle responsible for inflation, fluctuates with a 
minimal amount of self interactions. In fact, such a small 
amount of self interactions ensures that the fluctuations 
are nearly Gaussian, and that any non-Gaussianity pro- 
duced would be too small for detection [2ll. [2^ . [23L [24. [25|] . 
Non-Gaussianity therefore would be a measure of either 
interactions of the inflaton [2(1 [2?| or any non-linearities 
[H, HI], an d a detection of non-Gaussianity would indi- 
cate a violation of slow-roll inflation. 

In this spirit, models of non slow- roll inflation or 
alternatives to inflation have been proposed to gener- 
ate large, measurable non-Gaussianities. The curvaton 
mechanism produces curvature perturbations associated 
with the fluctuations of a light scalar field whose energy 
density is zero 30]. The inhomogeneous reheating sce- 
nario can produce non-Gaussianity through modulated 
reheating during the reheating stage [3l|. Using mul- 
tiple inflaton fields that are allowed to interact, those 
interactions can be used to source non-Gaussianity [32| . 
Lastly, warm inflation [33[ , ghost inflation [34[ and string 
theory inspired D-cceleration [35| and Dirac-Born-Infeld 
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(DBI) inflation [36j models also give rise to a large non- 
Gaussianity (see review in Ref. [371]). 

To connect with observable measurements, the asso- 
ciated non-Gaussianity of the CMB can be described in 
terms of the second-order correction to the curvature per- 
turbations in position space with 

*(x) = L (x) + /nl |#(x) - (0l(x)) 2 ] , (1) 

where the non-Gaussianity parameter /nl describes the 
amplitude of the second-order correction. This form was 
first suggested by Salopek & Bond [H, HH to describe the 
non-Gaussianity in primordial perturbations from infla- 
tion and has been the subject of experimental constraints 
using a variety of CMB and large-scale structure data in 
recent years. 

Instead of constraints on the non-Gaussianity param- 
eter in the position space, recent studies make use of 
the bispectrum involving a three-point correlation func- 
tion in Fourier or multipole space. The configura- 
tion dependence of the bispectrum B{k\, &2, k%) with 
lengths (k\, tt2, k%) that form a triangle in Fourier space 
can be used to separate various mechanisms for non- 
Gaussianities, depending on the effectiveness of of the 
estimator used. To summarize the status of the non- 
Gaussianity measurements, an analysis with WMAP 3- 
year data first suggested a hint of a non-Gaussianity in 
the local model with 27 < / NL < 147 (95% CL), far 
above the value of /nl < 1 expected in simple, single 
field, slow-roll inflation models \3<^. The WMAP team's 
preferred measurement of non-Gaussianity parameter in 
5-year V and W-band data is -9 < /nl < HI (95% 
CL) pjj. The most recent constraint on /nl comes from 
studying the WMAP 5-year data with an optimal esti- 
mator leading to -4 < /nl < 80 (95% CL) [40j]. At the 
68% confidence level, with a value of /nl = 38 ±21, there 
is still some marginal evidence for a non-zero value of the 
non-Gaussianity parameter. If such a result were to con- 
tinue to hold with Planck, which increases the precision 
of /nl measurement by a factor of 3 to 4, then our simple 
inflationary picture would need to be revised to include 
a more complex model. 

In this paper, we will pursue a new measurement of 
the primordial non-Gaussianity parameter with a new 
estimator that preserves some angular dependence of 
the bispectrum. On the contrary, the estimators em- 
ployed by most CMB non-Gaussianity studies, including 
those by the WAMP team pjj], involves a measurement 
that compresses all information of the bispectrum to a 
single number called the cross-skewness computed with 
two weighted maps. Such a drastic compression lim- 
its the ability to study the angular dependence of the 
non-Gaussian signal and to separate any confusing fore- 
grounds from the primordial non-Gaussianity. In addi- 
tion to Galactic foregrounds, non-Gaussianity measure- 
ments could also be contaminated by unresolved point 
sources, mainly radio and dusty galaxies, and Sunyaev- 
Zel'dovich (SZ) clusters, among others [4l|. Given the 
increase in size of CMB data, especially with Planck, it 



is also necessary to develop accurate measurement tech- 
niques to extract /nl that are unbiased. 

Our estimator for non-Gaussianity uses a weighted 
version of the squared temperature-temperature angular 
power spectrum [42l [43| , which we refer to as the skew- 
ness power spectrum. This power spectrum extracts in- 
formation from the bispectrum as a function of the mul- 
tipole of one triangle length in the harmonic space, while 
summing all configurations given by the other two side 
lengths. The difference in spatial dependence based on 
how the maps are weighted provides ways to separate 
primordial non-Gaussianity from that of the foregrounds. 
Here, we account for both point source and lensing bis- 
pectra with latter resulting from the correlation of the 
lensing potential with secondary anisotropies [44l |45| . 

To summarize our main results, after marginalizing 
over the normalizations of point source and lensing- 
secondary bispectra, with the combination of V and W- 
band maps we are able to constrain /nl = 11.0 ± 23.7 
at the 68% confidence level or —36.4 < /nl < 58.4 at 
the 95% confidence level. We find that /nl is never 
incompatible with zero at 68% confidence when /nl is 
estimated in independent bins of width 200 between 
2 < I < 600. We find a significant contribution from un- 
resolved point sources, but failed to detect the lensing- 
secondary cross-correlations using the two statistics we 
considered here. 

In section £|TT] we review the background theory and 
in £11111 we review the estimator used and our simulation 
procedure to compute the uncertainties. In i ]IVI we dis- 
cuss our methods for analyzing and simulating data. In 
section <|V]we discuss our results. In section VII we con- 
clude with a summary of our results. 

II. THEORY 

To begin, we define multipole moments of the temper- 
ature map through 

a lm = J dhT(h)Yr*(n). (2) 

The angular power spectrum and bispectrum are defined 
in the usual way such that 

( a l imi a hm2) = ^i/roimj^i > (3) 
(ai imi ai 2 m 2 ai 3m3 } = f ^ ^ B hhh . (4) 

Here the quantity in parentheses is the Wigner-3j sym- 
bol. The orthonormality relation for Wigner-3j symbol 
implies 

Bllhh = ^2 ( mj m 2 ^ 3 J ( a 'i»H a i2m2 a i3m 3 ) -(5) 
m\m 2 m 3 ^ ' 

The angular bispectrum, Bi ± i 2 i a , contains all the informa- 
tion available from the three-point correlation function. 
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FIG. 1: Qf;(r) and Pi(r) with respect to / for r values defined as 
followed: r = c(to — ar) where to is the present day conformal 
time and cr = 235 Mpc. In these plots, a = 0.04, 0.4, 1.0 and 
1.8. 



For example, the skewness, the pseudo-collapsed three- 
point function of Ref. 14611 and the equilateral configura- 
tion statistic of Ref. [47| can all be expressed as linear 
combinations of the bispectrum terms (see Ref. [29| for 
explicit expressions and Ref. [48] for an expression relat- 
ing skewness in terms of the bispectrum). 



and r is the comoving radial coordinate. 
The two functions in are gi ven by 

Q; ( r ) = ~~ J k 2 dkg T i(k)ji(kr), 
A(r) = £ I ' tfdkPtWgriWMkr) . 



(8) 
(9) 



Here, -P$(fc) cx fc" s_4 is the primordial power spectrum 
of Bardeen's curvature perturbations, and gTi(k) is the 
radiation transfer function that gives the angular power 
spectrum as Ci = (2/tt) J k 2 dkP^,(k)g 2 -, l (k). In Fig. 1, we 
show four example cases of a(r) and /3(r). We generate 
them using a modified version of the CMBFAST code 
p9| and for our fiducial cosmological parameter values, 
consistent with WMAP 5-year best-fit model, as summa- 
rized in Table I. 



Parameter 


Value 


H 


71.9 km/s/Mpc 


tl b h 2 


0.02273 


n c h 2 


0.1099 


n s 


0.963 


r 


0.087 


A 2 


2.41 x 10" 9 


(Jo for Q 


2.197 mK 


an for V 


3.133 mK 


(j for W 


6.538 mK 


/sky 


0.718 



TABLE I: Cosmological and noise parameters used in our 
analysis. The first set is our fiducial cosmology model taken to 
be consistent with WMAP 5-year best-fit cosmology [12J . The 
second set of numbers is the normalization parameters related 
to the instrumental noise in each of the three frequency bands 
used for the analysis. / s k y is the fraction of sky unmasked by 
KQ75 mask. 



A. Primordial Non-Gaussianity 

Here we focus on the local form of the primordial non- 
Gaussianity. Using the second order correction to the 
curvature perturbations <!> in equation ([1]) and following 
the derivation in Ref. [38| . we write the angular bispec- 
trum of temperature anisotropics as 

POO 

= 2I hl2h / r 2 dr [a h (r)A 2 (r)/?„ (r) + (Perm.)] , 

JQ 

(6) 

where 

T = / (2/i + l)(2i 2 + l)(2/ 3 + l) (h h h\ 
hhh ~ V 4tt ^ j ' 

(7) 



B. Unresolved Point Sources 

In addition to the primordial bispectrum, we also ac- 
count for the non-Gaussianity generated by unresolved 
radio point sources. If the sources are Poisson dis- 
tributed, the bispectrum takes a simple from [!U with 

B hhh = ^i/ 2 /3 6 PS' ( 10 ) 

where 

b^^g^J^S^dS, (11) 

where dn/dS is the number counts of sources and g(x) 
maps flux density to thermodynamic temperature with 
g(x) = c 2 (e x - lf/2k B v 2 x 2 e x with x = his/k B T CMB w 
i//56.84GHz. This conversion can be simplified to g(x) = 
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liK/(99.273ysr- 1 ){e x - l) 2 /x i e x . When model fitting to 
data, we will ignore the exact number counts of the un- 
resolved sources and parameterize the uncertainty with 
an overall normalization 



x 1CT 25 sr 2 



(12) 



where the index i is for the three bands from WMAP (Q, 
V, and W) we use here. 

Here, we only account for the shot-noise contribu- 
tion from point sources, similar to the analysis of non- 
Gaussianity measurements by the WMAP team. It is 
likely that unresolved point sources are clustered on the 
sky, though existing WMAP data with measurements at 
the two-point function level only lead to an upper limit 
on the clustering amplitude of point sources [5fJ. In fu- 
ture, especially for non-Gaussianity measurement with 
Planck, it may be necessary to include the bispectrum 
generated by clustered point sources. 



C. CMB Lensing-Secondary Correlation 

The gravitational lensing effect of the CMB also gen- 
erates a bispectrum through correlations of the lensing 
potential with secondary anisotropies that are generated 
at late times 0, HH] . 

To understand this signal, we note that the lensed tem- 
perature fluctuation in a given direction is the sum of the 
primary fluctuation in a different direction plus the sec- 
ondary anisotropy 

T(n) = T p (n + V6) +T s (n) (13) 
« [^n + a L)Yi n (n) + af m 

lm 



xV6(n) • VFj m (n) 



Integrating by parts and simplifying further following 
leads to a bispectrum of the form: 



lens— sec 



Bf 




(2i 1 + l)(2I 2 + l)(2/ 3 + l) 



i 3 (i3 + 1) 



4tt 

+ Perm. 



(16) 



When calculating the CMB lensing potential- 
secondary anisotropy cross-correlation bf we will include 
both the integrated Sachs- Wolfe (ISW) and the Sunyaev- 
Zel'dovich (SZ) effects, with the latter modeled using 
the halo approach [5l|, H2, [Hj]. We will take the sum 
of the two effects such that bf = bj sw + bf z . The 
cross-correlation between lensin g p otential and ISW is 
calculated in the standard way [54[ l55j for the fiducial 
ACDM cosmological model, using only the linear theory 
potential. For the lcnsing-SZ correlation, the linear 
halo model takes into account the SZ profile obtained 
analytically in Ref. [56j combined with the halo mass 
function similar to calculations of the SZ angular power 
spectrum. When model fitting the data, we will param- 
eterize the overall uncertainty with a parameter rji for 
each of the WMAP bands such that B} e f s r sec oc m. 

While lensing modification to CMB bispectrum alone 
is not expected to make a significant correction to the 
non-Gaussianity measurement, analytical calculations of 
the lensing effect on the CMB bispectrum suggest that 
the lensing-secondary correlation will be the main con- 
tamination to a reliable measurement of the primordial 
non-Gaussianity parameter [4lll57ll58. 59]. This includes 
the lensing-ISW effect since SZ can be "cleaned out" in 
multi-frequency data such as those expected from Planck 
[48l l60l ]. It is due to this reason that we include the 
lensing-secondary correlation here. 



p s 

a lm + a lm 



I'm' 

x / dnY7 n *(n)Ve(n) • VY^'(n) . (14) 



ay, 



Utilizing the definition of the bispectrum in Eq. (JSJ), we 
obtain 



plens— sec 
n hhh 



E 



X 



h h h 
mi m,2 Trig 



drh / dnY™**(rh)Y^*(n)C h 



x 



VF™i*(m) • (Ve(rh)T s (n)) + Pcrm.(15) 



where the extra five permutations are with respect to the 
ordering of (h,l2,h). 



III. ESTIMATORS OF /nl 

We will now motivate a new estimator for measuring 
/nl- For this we introduce the squared temperature- 
temperature angular power spectrum and discuss its use 
as a probe of the angular bispectrum. We motivate a 
new estimator by revising the original form in Ref. 

Through the expansion of the temperature 



42] 



T(n) = 5> m y ; m (n), 



we can write 



2 

In 



J dnT 2 (n)Y7"*(n) 



(17) 



(18) 



We emphasize here that af m denotes the multipole mo- 
ments of the temperature squared map and not the 
square of the multipole moments of the temperature map. 
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FIG. 2: Contributions to expected from primordial non- 

Gaussianity and unresolved point sources. We show the case 
with /nl = 1 for primordial non-Gaussianity (top) , shot- noise 
from unresolved point sources with b pB — 1 (middle), and 
lensing-secondary signal with r/i — 1 (bottom). 



We can now construct the angular power spectrum of 
squared temperature and temperature as 



cf- 1 = 



21 



— V 

+ 1^ 



a lrn a lrn • 



(19) 



After some tedious, but straightforward algebra we can 
write the relation between the bispectrum of the temper- 
ature field and the angular power spectrum of squared 
temperature and temperature as 



— TB h 
+ 1 




(20) 



Here, we have made use of the relation 




mi m2 m 



h h V 
mi 77i2 m 1 



°U'°mm' 
21 + 1 



(21) 



FIG. 3: Contributions to -E ; 2_1 expected from primordial non- 
Gaussianity and unresolved point sources. We show the case 
with /nl = 1 for primordial non-Gaussianity (top) , shot-noise 
from unresolved point sources with bp B = 1 (middle), and 
lensing-secondary signal with r\i = 1 (bottom). Note the large 
difference in the y-axis scale from top curve Involving primor- 
dial non-Gaussianity to middle and bottom curves with point 
sources. As is known, the skewness power spectrum associ- 
ated with E maps is more sensitive to shot-noise bispectrum 
from point sources. 



As is clear C l ~ sums up all triangle configurations of 
the bispectrum at each of the side length I of the triangle 
in multipolar space. 

If a priori known that certain triangular configurations 
contribute to the bispectrum significantly one can com- 
pute this sum by appropriately weighting the multipolc 
coefficients. This is essentially what can be achieved with 
the introduction of an appropriate weight or a window 
function in equation (|18p . Though the analytical expres- 
sion for the two-to-one angular power spectrum involves 
a sum over the two sides of the angular bispectrum, the 
experimental measurement is straightforward: one con- 
struct the power spectrum by squaring the temperature 
field, in real space, and using the Fourier transforms of 
squared temperature values and the temperature field, 
with any weighting as necessary. 

This simple form of the skewness power spectrum has 
already been used by Szapudi & Chen (6l| to constrain 
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/nl = 22±52 (la) with WMAP 3-year data. The form of 
the skewness power spectrum as written exactly in equa- 
tion (|2ip is not useful for a primordial non-Gaussianity 
measurement. We describe how to filter data for a mea- 
surement of primordial non-Gaussinity below. 

A. Skewness Estimator 

To obtain a more useful form, it is useful to review the 
form of the skewness statistic employed by the WMAP 
team, which is originating from Ref. [62]. The skewness 
statistic makes use of two set of maps of the CMB sky as 
a function of the radial distance r: 

A(r,A) = J2 Y i™(n)A lrn (r) (22) 

Int 

B(r,A) = ^2Yi m (n)Bi m {r) , (23) 

I'm 

where 

A lm (r) = ^ha lm (24) 
0r(r) 

B lrn {r) = t^J. biaim . (25) 

Here Ci = Cibf+Ni where 6/ are the frequency dependent 
beam transfer functions and Ni is the power spectrum 
from associated simulated noise maps. We discuss both 
these quantities later. 

In A and B maps weights are such that they are are 
constructed from the theoretical CMB power spectrum 
Ci under the assumed cosmological model, the experi- 
mental beam bi, and the primordial non-Gaussianity pro- 
jection functions a/(r) and 0i(r) where ai(r) and 0i(r) 
are defined in equations (jHJl and ©. 

The WMAP team's estimator [62J uses an integration 
in the radial coordinate to obtain the skewness of the 
product of the A and B 2 maps 

SW = J r 2 dr J dhA{r,n)B 2 (r,h) . (26) 

In practice this skewness is corrected by an additional lin- 
ear term that corrects approximately the effects of partial 
sky coverage associated with the mask and non-uniform 
noise. This term is computed by combining observed map 
with simulated maps that are Monte-Carlo averaged (see 
Appendix A of Ref. 

As is clear from above S3 involves a complete com- 
pression of data to a single number. While in principle 
different sources of non-Gaussianities contribute to S3 
with a single number alone it is impossible to separate 
out the primordial value from the non-Gaussianities gen- 
erated by secondary anisotropics and other foregrounds. 
To some extent the separation is aided by a different set 
of maps that are weighted differently than the case of A 
and B maps. 
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10 100 200 300 400 500 600 

I 

FIG. 4: Expected error for /nl calculated based on the Fisher 
matrix approach for each of the three noise curves for the 
WMAP in Q, V, and W-bands and with / sky = 0.718 when 
using KQ75 mask. The Cramer-Rao bound ranges from about 
~ 21 in V-band to ~ 23 in Q-band. This estimate assumes 
that only the primordial non-Gaussianity signal is present in 
the bispectrum and ignores the degeneracies between primor- 
dial non-Gaussianity and other parameters, such as those re- 
lated to unresolved point sources. 

A map optimized for the non-Gaussianity of the form 
generated by shot-noise from point sources is the E map: 

E(n) = ^Y lm (n)Ei m (r) , (27) 

lm 

where 

E lm (r) = ^a lm . (28) 

Similar to Sab 2 > one can a ^ so compute a skewness associ- 
ated with E maps by taking S E 3 = J dnE 3 (h). WMAP 
team used the latter to constrain the normalization of 
the point source Poisson term with 6 ps . 



B. Revised Skewness Power Spectrum 

In order to revise the previously discussed skewness 
power spectrum, instead of simply integrating over the 
A and B 2 maps, we extract the multipole moments of 
the B 2 map and the product AB maps 

(B 2 ) lm (r) = J dnB 2 (r,n)Y lm (n) 

(AB) lm (r) = fd&A{r,A)B(r,&)Y lm (&). (29) 
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These two multipole moments then allow us to write the 
new skewness power spectrum appropriately weighted in 
the same manner as the previous skewness estimator: 



C?- 1 = (C?' B + 2C? B > B ) 



c, 



c, 



1 



; 

A,B 2 

21 ■ 1 

;;. \ it J_ 

i 



21 + 1 



" 2 dr 



r dr 



(30) 

^Real{A, m (r)(B 2 ) im (r)} 

m 

J2^nBim(r)(AB) lm (r)} 



To see how C 2_1 probes the primordial bispectrum, we 
can write the multipole moments of the squared B map 
as 



(31) 



I'm' fin" 



Pv(r)0i»(r) /(2I + l)(2I' + l)(2i« + l) 



4tt 



I V I" 




I V I" \, , 

I „ a l'm' a l"m" J 
777 777 777 / 



where a[ m are the beam times the observed multipole 
moments (hai m ). Note that the observed multipole mo- 
ments relate to theory moments via another beam factor. 

Similarly, the multipole moments of the (AB) product 
map is 



{AB) lm (r) 



(32) 



EE 



ai>(r)0i»(r) /(2i + l){2l> + 1)(2Z" + 1) 
Cv C v , 



4tt 



I'm' l"m'' 

( I V l"\ ( I V I" 
1 J I m rri m" / 

The C l ' power spectrum is simply then 



m I'm' l"m" 



I3 v (r) f3i„(r) ai(r) (21 + 1)(2Z' + 1){21" + 1) 



C v C v , Ci V In 

( I V I" \ ( I V I" 

\ J \ 777 777' 777 /; 

Using the definition of the angular bispectrum, we can 
simplify to obtain 



21- 



- f r-dr^ 
J VI" 



(34) 



Cv Cv Ci 
' I V I" ^ - 



v {r) /3,«(r) a,(r) (21 + l)(2l> + l)(2l" + 1) 



4tt 







B'ivv'bib\b'l , 



(35) 



where B'nnn is the bispectrum estimated from data un- 
der beam smoothing. It relates to the theory bispectrum 
Bu> i" as B[ vv , — Biviiibib[b" . 

We can similarly simplify the term for Cf B ' B and 
putting the two terms together, we find that the total 
is 



C 



(3 



A.B* 



2C AB, B) 



(36) 



(21 + 1) 



E 



B, 



JVG,/nl = 1 
U'l" 



B'uniMm 



CiCi'Cv 



If we assume that the observed bispectrum is sim- 
ply that of the primordial non-Gaussianity then Biw = 
/nlSS/w and we can write an estimator for /nl as 



f NL = (21 + l)Cr 7Fng,ng(0 , (37) 

where .Fng,ng(0 is simply the Fisher matrix element for 
the primordial bispectrum with /nl = 1: 



Bivi"B\ vv , 
C'iC'vC'v- 



(38) 



where now we have redefined noise to be such that C'i 
( ': ■ X: ii as the bispectra are no longer beam smoothed. 

In reality C ; 2_1 includes contributions from secondary 
anisotropics and foregrounds. Here, we include the non- 
Gaussianities generated by point sources and the lensing- 
secondary correlation. Thus, we write 

(2l + \)Cf- 1 = /iVLF N G,NG(0 + ^ NG> ps + r,F NG 

,lcn— sec j 

and consider a joint estimation of the three unknown 
parameters. 

To help break degeneracies between the three param- 
eters, we also estimate the skewness power spectrum of 
the E map defined in equation P7| as 



C, 



E.E 



1 



2Z + 1 



^Real{£; ; , 



(40) 



(33) Similar to our derivation above one can simplify the mul- 
tipole moments of the (E 2 )i m to show that this probes 



E, 



2-1 _ 



= c, 



(41) 



(21 + 1) 



E 

VI" 



B 



PS,b ps =l 
U'l" 



B'lwMb 



i 



CiCi'Cv 



Thus, we write 



(21 + l)Ef 1 = fNLFps,NG(l) + AFps.SN + f/Fps.lcns-scc • 

(42) 

The two equations (|39| and (|42|) will form the main 
set of equations that we will solve with our measure- 
ments. While we have not explicitly stated so far, these 
two quantities will be measured in 3 WMAP frequency 
channels making use of Q, V, and W-band data. We al- 
low for frequency dependence in A and 77, but assume 
/nl is the same independent of the frequency in all three 
channels. 
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FIG. 5: A maps for V frequency band. From upper left hand corner moving clockwise: r = 0.04, 0.4, 1.0, 1.8 



C. Approximate corrections for partial sky 



Similarly, for E l we find 



Before we move onto discuss data analysis and our sim- 
ulations to compute the covariances, we note that we also 
make a correction to both Cf~ l and Ef~ l to account for 
partial sky coverage and inhomogeneous noise. This is 
done in an approximate manner by making use of the 
equivalent form of the linear terms of the skewness statis- 
tic in the language of our skewness power spectrum. For 
the case of C 2_1 estimator the correction is derived in 
Ref. 0: 



_ 1 {c A ' b2 - 2C ( - A ' B ' )B - C A ^ B ^' 

fsh-i l 



sky 

AB.B _ r ,(AB),B _ £,B(A,B) _ n A(B,B) 



sky 



C? D >° - c, 



c 



>} (43) 

where f s k y is the sky fraction observed. The new terms 
are defined as, for example, 

i m 

(44) 

where i runs over a set of N simulations and {B D A )\ m 
are the coefficients of spherical harmonics for the map 
produced by multiplying the z th simulated A map with 
the B map derived from raw data. 



E ?- 1 = T-{ C i 

T.ski, { 



E,E 2 qE,(E 



2C, 



(E,E)E 



}, (45) 



where terms such as C E '^ E ^ can be written similar to 
equation (|44[) above with the replacement of E maps in- 
stead of A and B maps. 



D. Theoretical expectation 

In Figure[5]and Figure[3]we show the theoretical expec- 
tations for C, 2-1 and -E, 2-1 , respectively. We plot these 
for the Q, V and W band by making use of the beam func- 
tions bi and noise power spectrum estimate Ni that are 
described in Section IIVB 11 Here, we show the cases of 
primordial non-Gaussianity with /nl = 1, point sources 
with Ai = 1 and lensing-secondary cross-correlation with 
bf calculated for the sum of ISW and SZ effects with 
m = 1. 

As is clear from Fig. 2, the primordial non-Gaussianity 
signal is expected to be degenerate with foreground non- 
Gaussianities. The shape of Cf~ l alone is not enough to 
clearly separate primordial non-Gaussianity signal from 
point source and lensing non-Gaussianities. Fortunately, 
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FIG. 6: B maps for V frequency band. From upper left hand corner moving clockwise: r = 0.04, 0.4, 1.0, 1.8 



the separation is aided when C l is combined with 
Ef -1 . As is clear from Fig. 3 (espcially note the dif- 
ference in the y-axis range for the top and middle plots), 
this latter power spectrum allows a better determination 
of the point sources. In practice, we perform a com- 
bined analysis of both spectra, including the confusion 
from secondary bispectra, when model fitting to quanti- 
ties, /nl,^4i and r\i. To compare with previous results on 
the literature related to the non-Gaussianity parameter 
with WMAP data using the effects of point sources only, 
we also consider the case where lensing is ignored in the 
analysis. 

In Fig. 4, we include a plot of the ^ng ng(0; show- 
ing the expected error /nl as a function of the multi- 
pole for Q, V, and W bands. Out to l max of 600 and 
with / s ky = 0.718, the Carmer-Rao bound is at the 
level of ~ 21 with V-band to ~ 23 with Q-band. This 
assumes that bispectrum only contains primordial non- 
Gaussianity, but the degeneracy between secondary non- 
Gaussian signals and primordial non-Gaussianity is ex- 
pected to increase the optimal error at some level more 
than this bound. Also, to saturate the Cramer-Rao 
bound an optimal estimator that accounts for the mode- 
mode correlations associated with the partial sky and the 
mask will become necessary [40j . Our estimator only ac- 
counts for the cut-sky approximately making use of the 



linear terms. We also weight each multipole coefficient 
with (Citif + Ni) -1 , as in the case of Gaussian statistics 
appropriate for the whole sky. While this approach is 
not different from that of the WMAP team's [l2| , in an 
upcoming paper we hope to return to the issue of an ex- 
act calculation implementing the full covariance for the 
two-to-one skewness power spectrum. 



IV. DATA ANALYSIS 

We first discuss our data analysis procedure and then 
how we computed the covariance through simulations. 



A. Measurement of Cf 1 and E? 



To extract Cf 



_1 and Ef~ from data we use the raw 
WMAP 5-year Stokes-I sky maps for the Q, V and W fre- 
quency bands as available from the public lambda web- 
site 1 . We use Healpix 2 [63j to analyze the maps. Specif- 



1 http://lambda.gfsc.nasa.gov 

2 For more information sec http://hcalpix.jpl.nasa.gov 
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FIG. 7: E maps for Q, V and W frequencies. 



ically, starting from the fits files of raw maps we use 
anafast, masking with the KQ75 mask and without an it- 
eration scheme, to generate multipole coefficients (a; m s) 
for each frequency map out to Z max = 600. We will refer 
to these multipole moments hereafter as as af m . With 
these definitions we use equations (f!M]) . (j2"5")l , and (|28p to 
generate A\ mi Bi m , and E\ m by substituting af m in place 
of a lm . 

Our recipe for obtaining U?~~ and Ef' 1 is: 

1. Use Healpix and the KQ75 mask to generate af m 
from the WMAP 5-year Stokes-I Sky Maps for the 
Q, V and W frequency bands. 

2. Obtain A(r,h), B(r,h) and E(h) from equa- 
tions (Jill), (gU, and (US]) by using a(r) and (3{r) 
from equation © and ©, respectively, with the 
replacement of aim with af m . Figures O M an d [7] 
show the resulting maps. 

3. Calculate Cf' B , C B ' AB and the linear terms from 
equations ([ST))) and (|4"3")l . respectively, with latter 
using equations of the form (|4"4")) . Repeat the same 

to obtain C E ' E with E maps as defined in equa- 
tion (140|) and the corresponding equations for linear 
terms in equation (145|) . These correction associated 



with partial sky coverage involves the use of sim- 
ulated maps described below. We integrate over r 
from t = 0.004 to 2 with 500 steps, (see Fig. 1). 

4. Use the estimated Cf " 1 and Ef" 1 with WMAP Q, 
V, and W maps for our parameter estimate analysis 
(see below). 

5. Compute analytically Fij terms with / max = 600 in 
each of the summations of l\ and li and making 
use of the noise and beam spectra for WMAP (see 
below) . 

In Figure [5] we see the Cf^ 1 for each WMAP frequency 
band plotted as a function of /. These plots were gener- 
ated by binning the estimators with 51 of 40 and plotting 
the midpoint of each bin. The V and W estimators have 
roughly the same shape and are mostly positive. The Q 
estimator is noticeably different, dropping negative when 
I > 300. 

Furthermore, in Figure [5] we see Ef- 1 for each WMAP 
frequency band plotted as a function of /. Like the esti- 
mators mentioned above, these were similarly binned in 
bins of size 61 = 40. 

Lastly, in Figures [T3] and [TJ] we see the contributions 
to Cf 1 and Ef -1 from each term in equations H31 and 
US] respectively. The linear terms are not very significant 
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FIG. 8: Of' 1 (top) Ef- 1 (bottom) for Q, V and W with the 
measured spectra binned with a width of 81 = 40. 
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FIG. 9: The noise power spectrum from our masked simulated 
noise maps for the Q, V and W frequency bands compared 
with the analytical results from the WMAP 5-year team. 




FIG. 10: The beam transfer functions, bi, used in our analysis 
for each frequency Q, V and W. 



compared to the other terms. Nevertheless, they are still 
considered in this analysis. 



B. Simulations of A, B, and E maps 

In order to do proper statistics for Cf~ l and Ef , 
we create 250 simulated maps at each WMAP frequency 
band. To do so we first produce 250 Gaussian maps to 
model the CMB sky. For the Gaussian maps we run 
synfast routine of Healpix with an in-file representing the 
WMAP 5-year best-fit CMB anisotropy power spectrum 
and generate maps with information out to I — 600. We 
then use anafast, without employing an iteration scheme, 
masking with the KQ75 mask, to produce a/ m 's for the 
Gaussian maps out to I = 600. We will refer to these 
aim's now collectively as af m . 



1. Noise 

In addition to these Gaussian maps we create 250 noise 
maps for each each of our frequency bands: Q, V and W. 
We generate these maps from white noise with mean = 
and standard deviation = 1 taking into account oo & n d 
Nobs as follows: 



N(h) 



(TO 



zn(h) 



(46) 



obs 



where N(h) is our noise map and n(n) is a map made 
of pure white noise, A^bs is the number of observations 
per pixel and <tq is the rms noise per observation. We 
use the frequency dependent A^bs for each point in the 
sky provided by the WMAP 5-year Stokes-I map fits files 
and take <tq = 2.197, 3.133, and 6.538 rnK as established 
by the WMAP team for Q, V and W band 5-year data 
respectively [gl, [6^| . See also Table HI 
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FIG. 11: The results for Cf -1 from all 250 simulations for the 
three frequency bands. These plots have been binned with 
SI = 10 



Starting with these noise maps we create alms using 
anafast with the KQ75 mask with no iteration scheme 
out to I = 600. We will henceforth refer to these alms 
collectively as afL. Furthermore, to calculate the power 
spectrum from these noise maps we use Healpix to eval- 
uate the analytical expression: 



N t =0 



pix 



d 2 n a$M(n) 
4vr/ sky AUs(n) ' 



(47) 



where J7pi x = 47r/iV p i X is the solid angle per pixel, M(n) 
is the KQ75 mask, / s k y = 0.718 is the fraction of sky 
retained by the KQ75 maskfl^. Fig. [9] shows the power 
spectrum from our simulated noise maps for each fre- 
quency compared to the analytical values quoted by the 
WMAP 5-year team [12]. For reference, the beam func- 
tions b\ are plotted is Fig. [I0l 

To use our estimator on the simulated maps we must 
add the noise to the Gaussian maps while at the same 
time correcting for the beam. To do this we work in 
multipole space where we construct the total simulated 
af m = af m bi + a^ m where bi are the frequency dependent 
beam transfer functions plotted in Fig. fTTJl 

FigurefTTlshow the results of C 2 " 1 plotted with respect 
to I for each frequency band. Similarly, Figure [T^] show 




FIG. 12: The results for £f -1 fxom all 250 simulations for the 
three frequency bands. These plots have been binned with 
61 = 10 



the all 250 simulated -Ef -1 plotted for each frequency 
band. These were binned with SI = 10. 

From these 250 simulations we are able to develop a 
covariance matrix that will be used for best fit estimates 
with error bars. We find this covariance matrix by bin- 
ning all 250 resulting estimators, Cf" 1 (or Ef" 1 ), in bins 
of 51 — 40. We can then treat each of these as an obser- 
vation for each I bin and create the covariance matrix 
by calculating the covariance of these observations. This 
produced and NxN covariance matrix where N is the 
number of I bins. 

Figure IT51 shows the correlation matrices from the sim- 
ulations. These matrices were obtained by taking the 
covariance matrix, Cy and building the correlation ma- 
trix Ca from the normalization: 



Cij — 



\JCaC 



(48) 



We see that the correlation matrix obtained from the 
Cf" 1 simulations show that these estimators have highly 
correlated I bins. It is interesting to note that the low I 
bins are highly correlated with each other and the high I 
bins are highly correlated with each other but low I bins 
are not correlated strongly with the high ell bins. 

We also see correlation in the Ef estimators but not 
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FIG. 13: Contributions to C\ 2 for Q, V and W maps 

nearly to as great a degree as C 2 ~ l above. Furthermore, 
in the bottom of Figure [18] we see the full Cf- 1 + E 2 ~ Y 
correlation matrix and note there is correlation between 
the I bins between C 2 ~ l and E 2 ~ l but not as much as 
there is between C, 2_1 alone. 



C. Best Fit Estimation 

In order to fit the data, we use a least squares fitting 
analysis. Given a data set consisting of n points (xi,yi), 
we can fit this data with a model function /(x, p) where 
there are m adjustable parameters held in the vector p. 
We wish to find which of those parameter values best fit 
the data. 

To do this we minimize the y 2 value defined as: 

X 2 = (y T -M-p) T C- 1 (y-M-p). (49) 

where y defines our data points we would like to fit to, 
p are the parameters we wish to solve for, M is a matrix 
containing our theoretical model we use for fitting and C 
is our covariance matrix described above. 

For example, for a single frequency analysis where we 
would like to fit for /nl and the coefficients for point 
sources: y = Cf~ 1D taken from data, M is the vector 
containing (Cf- 1Th , PS Th ) and p = (/ NL ,^> with A, 
being the coefficient for point sources. 
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FIG. 14: Contributions to E}~ 2 for Q, V and W maps 

We minimize \ 2 by setting its derivative to zero and 
solving for p yielding: 

p = (M T C" 1 M)" 1 M T C^ 1 • y. (50) 

Lastly, we find the error bars for our best fit parameters 
via 

Ap 2 = (MC^M)" 1 (51) 

where the diagonal of this matrix gives the variance of 
the parameters and the x 2 fit is given by equation (|49|) . 

V. RESULTS AND DISCUSSION 

A. /nl estimate 

We now discuss the results of our analysis. The pri- 
mordial and foreground non-Gaussianity parameter es- 
timates are summarized in Table |TT] for the case with 
and without point sources and in Table IIIII for the case 
with both point sources and lensing-secondary correla- 
tion. For each of these analyses, we bin and tabulate our 
measurements with bins of SI = 40. When SI < 20 the 
data are noisy to see the overall structure with a large 
covariance between adjacent bins and when SI > 100 in- 
formation from the fluctuating point source curves is lost 
leading to a large degeneracy between parameters and an 
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FIG. 15: The various /nl values taken from various binnings 
of the data from a full analysis considering both Cf^ 1 and 
Ef^ 1 with point sources and ISW. 
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FIG. 17: Histogram of the /nl estimated from Gaussian and 
noise simulations for the cross-skewness statistic S3 = "^2(21 + 
l)Cf _1 . Top: Q, Mid: V, Bottom: W. A best fit Gaussian 
curve is plotted in blue over each histogram. 
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FIG. 16: Angular dependance of /nl between 2 < I < 600 
in bins of 200. Top is for Cf' 1 only and with point sources. 
Middle is the full measurement with Cf~ x and -Ef -1 and using 
point sources. The lower panel is full measurement with both 
point sources and lensing-secondary correlations. The blue 
circles use mutually disjoint bins from each other. The green 
triangles also use mutually disjoint bins. 



increase in parameter errors. Furthermore, of all binning 
widths between 20 < SI < 100 the results are similar, 
but the best % 2 value is always found with a binning at 
SI = 40 (Fig. [TSJ. Note that in the limit of a large SI 
bin (with SI > 200) we effectively reach the case of de- 
termining /nl similar to the previous skewness statistic, 
with effectively just one data point per band. 

In Table II and III the first set of results, denoted by 
C?~ , show the case when we fit our measured Cf^ 1 
to the theoretical predictions involving a combination 
of primordial non-Gaussianity, point sources, and lens- 
ing correlations as shown in Fig. Given that Q 
map leads to a poor % 2 when model fitting Q alone 
or Q in combination with other maps, we exclude the 
Q+V+W combination and instead consider V+W as 
our preferred set of maps. When fitting to V and W, 
we compute the covariance of V and W, for example 



•,2-1, W\ 



^-X,V C 2-X,W ) _ {c 2-l,V ){c ,-i,vvy Without pomt 

sources and lensing and simply fitting to /nl with Cf^ 1 
we find 4.8 ± 27. 7 3 . If the shot-noise from point sources 



3 We quote la results with ± error and 2<r result as a range. 
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are included, after marginalizing over Ay and Aw, we 
find / NL = 39.0 ± 30.7. 
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FIG. 18: Matrices showing correlations between each fre- 
quency for both Cf -1 and E? -1 . The upper left hand corner 
shows correlations between small ell, and moving toward the 
lower right corner shows correlations between high ell. 

As we discussed earlier, however, fitting to C, alone 
with point sources lead to a worse determination of /nl 
than the case where point sources are ignored due to 
the degeneracy between primordial non-Gaussianity and 



point sources. Thus, we also include E maps in our anal- 
ysis with the associated results from the skewness power 
spectrum denoted with E 2 ^ 1 in Table |TT] and IIIII The 
E maps provide a better estimator for the point sources 
but a worse estimator for /nl than the C 2_1 estimator, 
for reasons we discussed already. (Fig. [3]) With E 2_1 , 
the error bars for the point source amplitudes are about 
half of what they were for C 2_1 alone, whereas the error 
bars on /nl are about three times worse. 

One interesting thing to note is that E 2_1 is always 
positive. This shows up in the best fit /nl values, where 
for Cf _1 only the Q map pushes /nl towards a negative 
value, whereas Ef~ l from Q pushes /nl to a large pos- 
itive value. In fact, if we include Q band and do a /nl 
analysis with the E map alone, we find a 6a detection 
of the primordial non-Gaussianity. The % 2 from such an 
analysis, however, is poor and the result should not be 
trusted as a detection of a non-zero /nl- 

Finally we consider the best fit when C 2_1 and E 2_1 
are combined. The V+W analysis gives us the best 
constraint on /nl with —36.4 < /nl < 58.4 at the 
95% confidence level or (11.0 ± 23.7) at the 68% con- 
fidence level, when we include both point sources and 
the lcnsing-secondary correlation and marginalize over 
(Ay , Aw , T]v , Vw) ■ As with the C ; 2-1 only analysis, this 
combined analysis has /nl consistent with zero at la. As 
can be seen by comparing Table II and III V+W case, 
our /nl is essentially the same whether we include the 
lensing-secondary bispectrum or not. 

While we do not include Q-band in our /nl estimate, 
the Q-band point source amplitude of (24.2±5.0) x 10~ 25 
sr 2 using the combination of C 2_1 and E 2_1 is consis- 
tent with the WMAP team's preferred value for t he p oint 
source amplitude of (4.3 ± 1.3) x 10~ 5 /iK 3 -sr 2 pj. In 
their units, our 6p S is equivalent to (4.9 ± 1.0) x 10~ D 
/iK 3 -sr 2 . While we cannot make an exact compari- 
son as WMAP team tabulates their point source val- 
ues with Z max of 900, our values for 6p S and b^ s are 
also within uncertainties consistent with previous mea- 
surements. While the non-Gaussianity associated with 
point sources is detected, we do not detect the lensing- 
secondary bispectrum. It is likely that the C?~ and 
E 2_1 are not the best ways to detect this correlation. 
The best-fit values for r\y and r\w, however, are close to 
their la errors. 

As tabulated in Table III, including the bispectrum of 
lensing-secondary correlations does not lead to a signif- 
icant degradation of /nl measurement. We find /nl = 
11.0 ±23.7 at the la confidence level, but we do not find 
a detection of r\i in each of the three bands when C 2_1 
is combined with E\ r . 

Note that our V+W analysis gives a /nl value fully 
consistent with zero at the la level. Previous results 
have suggested a marginal hint of a primordial non- 
Gaussianity with the most recent optimal anlaysis giving 
/ NL = 38+21 [H (Table V). Compared to this result, our 
V+W has a slightly worse error with 11.0 + 23.7, and the 
increase of 13% is consistent with the fact that our anal- 
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ysis is sub-optimal. As discussed earlier, however, our 
approach is not different from both the WMAP team's 
approach [l2j and previous other estimates of /nl [39l |- 
Moreover, our Z max is set at 600, while their analysis ex- 
tends to 750. 

To see if there is any scale dependence to non- 
Gaussianity we bin /nl hi widths of 200 and estimate 
the value between < I < 600. The results are shown 
in Fig. 16 and tabulated in Table IV. Except in the 
last bin for the case with point sources only between 
400 < I < 600, our / NL values are fully consistent with 
zero at the la level and the last bin is consistent with 
zero at the 2a level. The last bin also has a large error 
due to the increase of the instrumental noise. For the 
same reason, we do not pursue a measurement of /nl 
when Z > 600. 

It is also interesting to note how accurate our overall 
error estimate is. As we compute our covariances with 
250 simulations there is an inherent error of 1/V250 in 
the error bars we obtained in this analysis. Because of 
this, we note that a more accurate estimate of /nl should 
be to consider it as 11.0±23.7(±1.5) where the extra error 
within the bracket denotes an additional statistical error 
associated with the finite number of simulations. 



B. Cross-Skewness 

Previous results for /nl from the WMAP 5-year team 
compute /nl by compressing all information into a sin- 
gle quantity called cross-skewness defined by equationl26"l 
To compare our measurement C 2_1 with their results we 
calculate our own equivalent version of this cross skew- 
ness statistic defined as 

Sab- = £>i + l ) C i ~ 1,D ( 52 ) 

where Cf~ l ' D is the estimator obtained from data. We 
also compute the skewness of the E map using E?~^' D in 
above. We jointly fit Sab 2 and Se 3 with a combination 
of /nl and Ai by effectively comparing the statistic from 
data to prediction from theory with theory expectation 
computed as, for example, Sab 2 — + l)Cf~ 1 ' Th . In 
order to determine the errors we also preform the same 
cross-skewness analysis on all 250 simulations and cal- 
culate the covariance of Sab 2 an d Se$ from these 250 
numbers for each frequency. We find that /nl estimated 
from each of the 250 Gaussian and noise simulations lead 
to a Gaussian error distribution (Figure [T7| . 

We tabulate our results for /nl after marginalizing 
over AiS in Table VI. Here, when doing the summa- 
tions we set / max = 500 so we can compare directly with 
WMAP 5-year published results [HJ . We see that for all 
three channels we have good agreement with the WMAP 



team's 5-year findings. Our best-fit value tends to be 
bit more positive than quoted by the WMAP team (with 
0.26(7, 0.25ct, 0.17cr in Q, V, and W respectively), but this 
is a small difference when compared to the large error 
bar. The errors quoted in the WMAP 5-year paper is 
consistent with our measurements had we used the skew- 
ness statistic. However, as we discussed earlier, fitting 
to C 12-1 and E^~ x leads to an improvement in the error 
estimate of /nl since the shapes of the two skew spec- 
tra allow us to break the degeneracies better. Comparing 
our V+W result using the two spectra to skewness for the 
same maps, we find that the improvement in the error is 
roughly 20%. 



VI. CONCLUSION 

In this paper, we constrained the primordial non- 
Gaussianity parameter of the local model /nl using the 
skewness power spectrum associated with the two-to- 
one cumulant correlator of cosmic microwave background 
temperature anisotropies. This bispectrum-related skew- 
ness power spectrum was constructed after weighting the 
temperature maps with the appropriate window func- 
tions to form an estimator that probes the multipolar 
dependence of the underlying bispectrum associated with 
primordial non-Gaussianity. 

We also estimate a separate skewness power spectrum 
more sensitive to unresolved point sources. When com- 
pared to previous attempts at measuring the primordial 
non-Gaussianity with WMAP data, our estimators have 
the main advantage that we do not collapse information 
to a single number. When model fitting two-to-one skew- 
ness power spectrum we make use of bispectra generated 
by primordial non-Gaussianity, radio point sources, and 
lensing-secondary correlations. W 

We analyze Q, V and W-band WMAP 5-year data us- 
ing the KQ75 mask out to Z max = 600. Using V and 
W-band data and marginalizing over model parameters 
related to point sources, our overall and preferred con- 
straint on /nl is 11.0 ±23.7 at the 68% confidence level 
(—36.4 < /nl < 58.4 at 95% confidence). Despite pre- 
vious claims, we find no evidence for a non-zero value of 
/nl even marginally at the la level. 
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Type 


/nl (no PSs) 


/nl (w/PSs) 


Aq 


A v 


A w 


X 2 /dof 
















Q 


-61.6 + 32.2 


-10.5 + 33.6 


62.0 + 12.1 






1.6 


V 


5.4 + 30.4 


36.5 ± 32.9 




22.1 + 9.0 




0.6 


w 


5.5 + 30.0 


31.8 + 33.3 






18.6 ± 10.1 


0.6 


v+w 


4.8 + 27.7 


39.0 ± 30.7 




18.5 + 8.2 


25.3 ± 9.2 


1.0 


Ef- 1 














Q 


426.4 ± 100.5 


191.5 ± 115.4 


57.0 + 13.8 






1.3 


V 


159.1 + 98.0 


94.2 + 106.6 




13.8 + 8.9 




0.3 


w 


90.4+ 102.6 


49.2 + 112.4 






9.6 + 10.7 


0.3 


v+w 


133.1 ± 140.9 


69.8 + 100.6 




16.2 + 8.6 


9.4 + 10.3 


0.8 


Full 














Q 


-23.1 + 29.4 


-22.0 ± 29.4 


24.2 + 5.0 






3.2 


V 


13.1 + 26.8 


16.3 + 26.8 




4.2 ± 2.1 




0.6 


w 


19.5 + 26.9 


19.5 + 27.0 






0.4 + 2.3 


0.6 


v+w 


11.4 + 23.6 


11.5 ± 23.6 




5.0 + 1.8 


-1.8 + 2.0 


0.9 



TABLE II: Parameter estimates with Cf~ x (top), Ef' 1 (middle), and the combination of the two (bottom) with Q, V, W and 
V+W maps for the case where we ignore point sources and including point sources. The point source amplitudes are listed 
under columns for Ai's. The Q-band point source amplitude of (24.2 + 5.0) x 10~ 25 sr 2 , equivalent to (4.9 + 1.0) x 10~ 5 ^K 3 -sr 2 
is consistent with the WMAP team's preferred value of (4.3 ± 1.3) x 10~ 5 /iK 3 -sr 2 . The value of /nl with the amplitude of 
point sources marginalized over 11.5 ± 23.6. 



Type 


/nl (PS + lensing) 


Aq 


A v 


A w 


VQ 






X 2 /dof 




















Q 


21.1 + 40.3 


-80.2 + 39.3 






-11.7 + 5.8 






3.4 


V 


15.7 + 38.9 




8.7 + 23.0 






-3.7 + 4.6 




1.0 


w 


-13.5 + 39.8 






39.7 + 25.6 






0.6 + 4.4 


1.2 


v+w 


14.3 ± 37.6 




18.2 ± 20.8 


9.0 + 22.0 




-2.7 + 4.1 


-2.2 + 4.0 


1.3 


Ef- 1 


















Q 


122.2 ± 118.6 


8.5 + 6.2 






6.6 + 1.7 






0.7 


V 


80.5 ± 107.8 




2.1 ± 2.6 






1.2 ± 1.1 




0.3 


w 


62.3 ± 113.2 






-0.2 ± 2.5 






0.9 + 1.3 


0.3 


v+w 


72.0 ± 103.1 




1.9 + 2.4 


-0.5 + 2.4 




1.4 ± 1.1 


1.3 + 1.2 


0.8 


Full 


















Q 


21.8 + 29.6 


24.0 + 5.7 






0.2 + 1.2 






3.3 


V 


16.7 + 27.1 




4.1 ± 2.4 






0.2 + 0.5 




0.6 


w 


18.7 + 27.2 






0.5 + 2.3 






-0.3+1.0 


0.8 


v+w 


11.0 ± 23.7 




2.8 ± 2.2 


-0.4 + 2.2 




1.0 + 0.8 


-0.6 + 0.9 


0.9 



TABLE III: Parameter estimates with Cf' 1 (top), Ef -1 (middle), and the combination of the two (bottom) with Q, V, W 
and V+W maps for the case where we account for both point sources and the amplitude of lensing- secondary bispectrum. The 
point source amplitudes are listed under columns for Ai's, while the amplitude of lensing-secondary signal is tabulated under 
r\i 's. Our preferred value of /nl with the amplitude of point sources and the lensing-secondary signal marginalized over using 
V and W maps in combination is 11.0 ± 23.7. 
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